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ABSTRACT 

High angular resolution images obtained with a hypertelescope can strongly constrain the radiative-hydrodynamics 
simulations of red supergiant (RSG) stars, in terms of intensity contrast, granulation size and temporal variations of the 
convective motions that are visible on their surface. The characterization of the convective pattern in RSGs is crucial to 
solve the mass-loss mechanism which contributes heavily to the chemical enrichment of the Galaxy. We show here how 
the astrophysical objectives and the array configuration are highly dependent to design a hypertelescope. For a given field 
of view and a given resolution, there is a trade-off between the array geometry and the number of required telescopes 
to optimize either the (u,v) coverage (to recover the intensity distribution) or the dynamic range (to recover the intensity 
contrast). To obtain direct snapshot images of Betelgeuse with a hypertelescope, a regular and uniform layout of telescopes 
is the best array configuration to recover the intensity contrast and the distribution of both large and small granulation 
cells, but it requires a huge number of telescopes (several hundreds or thousands). An annular configuration allows a 
reasonable number of telescopes (lower than one hundred) to recover the spatial structures but it provides a low-contrast 
image. Concerning the design of a pupil densifier to combine all the beams, the photometric fluctuations are not critical 
(A photometry < 50%) contrary to the residual piston requirements (OPD < /1/8) which requires the development of an 
efficient cophasing system to fully exploit the imaging capability of a hypertelecope. 

Keywords: Interferometry, imaging, direct imaging, snapshot, hypertelescope, densification, stellar atmospheres, red su- 
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1. INTRODUCTION 

Current ground-based optical interferometers are already able to provide images, using the model fitting and the aper- 
ture synthesis techniques.-' 5 However, these images are limited both in spatial resolution and in spatial frequency coverage, 
so that the image reconstruction is always an undetermined problem and requires a pragmatic method to be achieved. Taking 
into account the astrophysical objectives in high angular resolution,- the future interferometric arrays should have many 
telescopes with kilometric baselines to enhance their imaging capabilities. They also require a suitable beam-combiner to 
manage the high number of entrance sub-apertures. Current pair-wise beam-combiner are not suitable in that case, and an 
all-in-one combination scheme should be used. 

The hypertelescope concept appears as a convenient and efficient solution.EH Using many telescopes, an efficient co- 
phasing system and a densified pupil combiner, they can provide direct snaphot images of complex targets. Due to the 
highly diluted entrance pupil of the telescopes spread on several hundred meters or kilometers, a multi-axial combination 
leads to a faint and spread Fizeau image formed in the common focal plane. A pupil densifier rearrange the beams so as 
to preserve the high resolution information and to highly intensify the signal regarding to the Fizeau mode. By increasing 
the relative size of the beams and keeping the general layout of the position of the beams, most of the light is concentrated 
in the useful field of view (FOV),JA reducing the image width to a convenient size for the detectors whatever the lengths 
of the baselines. The characteristics of the point spread function (PSF) have already been described in a previous paperpS 
showing that there is a trade-off between the angular resolution, the useful field of view and the halo level of the PSF to 
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choose an optimal array configuration. This paper deals with the design of an array and of a combiner in the case of stellar 
surface imaging applications. 

We concentrate here on the case of red supergiant (RSG) stars. Massive stars with masses between roughly 10 and 
25 M Q spend some time as RSGs being the largest stars in the universe. They have effective temperature (r e ff) ranging 
from 3450 to 4100 K), luminosities of 20000 to 300000 L , and radii up to 1500 R Q ^ which makes them some of the 
brightest stars known. Such extreme properties foretell the demise of a short-lived stellar king because they are nearing 
the end of their life and they are doomed to explode as a supernova. RSGs still hold several unsolved mysteries : (i) the 
mass-loss mechanism, shedding tremendous quantities of gas, is unidentified ; (ii) their chemical composition is largely 
unknown because of difficulties in analyzing their complex spectra due to the low surface temperatures and vigorous 
convection. RSGs contribute extensively to the chemical enrichment of our Galaxy loosing enormous quantities of their 
mass due to an unknown process. The vigorous convection that they experience could be at the base of the mass-loss 
mechanism. 

Recently, 3D radiative-hydrodynamics (RHD) simulatio ns don e with the numerical code C0 5 BOLD^^ coupled with 
the calculations done with radiative transfer code OPTIM3D^U2l have reported robust interferometric comparisons of hy- 
drodynamical simulations with existing observations^^ in the optical, H and K band regions, arguing for the presence of 
convective cells of various sizes on the red supergiant Betelgeuse. The characterization of the convective pattern in RSGs 
is crucial to solve the mass-loss problem and high angular resolution imaging with hypertelescopes will constrain the RHD 
simulations in terms of surface intensity contrast, granulation size and temporal variations. 



2. ASTROPHYSICAL REQUIREMENTS IN THE CASE OF BETELGEUSE 
2.1 Astrophysical parameters used for radiative-hydrodynamics (RHD) simulations 




Figure 1. Original map of the linear intensity of Betelgeuse at 7450 A, with 2 different scale to underline the structures, as the intensity 
covers a large dynamic range. The range is [0,600 000] erg/s/cm 2 /A(left) and [0,300 000] erg/s/cm 2 /A(middle), whereas the maximum 
intensity of the central bright cell reaches 1200 000 erg/s/cm 2 /A. Synthetic spectrum of Betelgeuse on the range [7100A-7800A] (right). 

Actual RHD simulations aim to reproduce the stellar parameters of the prototypical RSG star Betelgeuse (HD 39801, 
Ml-2Ia-Ibe), which is one of the best studied RSGs in term of multi-wavelength imaging because of its large luminosity 
and angular diameter. The model used in this work has the closest stellar parameters corresponding to this star : a 1 2 M 
stellar mass, a numerical resolution of 235 3 grid points with a step of 8.6 R Q , a luminosity of L = 93 000+ 1300L Q averaged 
over the spherical shells and over time, an effective temperature of T e s = 3490 + 13 K (compared to T e « = 3640 K EJ), a 
radius of R = 832 + 0.7 R Q , and a surface gravity log(g) = -0.337 ± 0.001 (compared to log(g) = -0.3 •% 

We have computed intensity maps using OPTIM3D for the RHD simulation at different wavelengths in the region bet- 
ween 710 and 780 nm (in total 700 images, 0.1 nm apart). This spectral region is interested by strong molecular absorption 
(namely TiO ; Fig. [D that cause high-contrast patterns characterized by dark spots and bright areas. The brightest areas 
exhibit an intensity 50 times brighter than the dark ones with strong changes over some weeksplas shown in Fig.Q] 



2.2 Astrophysical parameters for spectral analysis 

Three different spectral resolutions would be interesting to explore : 

- Low resolution (R tss 5) to study the large scale granule. This is more interesting in the H and K bands where only 
few convective cells of the size of « 60% of the stellar radius cover the whole surface of the simulated rsgP 

- Medium resolution (R a 20 - 30) to study the temperature stratification in the atmosphere. Observations at wave- 
lengths in a molecular band and in the close pseudo-continuum probe different atmospheric depths, and thus layers 
at different temperatures that highly constraint the limb-darkening effect of RHD simulations.^ 

- High resolution (R » 1000 or more) to study the cinematic in the spectral lines. In fact, these stars have velocities up 
to 20 - 30 km/s in the optical thin region that strongly affect the line broadening. 

In the following analysis of this paper, we consider only the intensity map at 745 nm (Fig.[T). 



2.3 Objectives to image Betelgeuse 

The main objectives that we need to recover from future interferometric images in order to better constrain the simula- 
tions are : 

1 . the size of the granulation cells from the large scale granule to the small-medium scale convective-related structures ; 

2. the surface intensity contrast ; 

3. the timescale of the convective motions. 

These objectives are crucial : 

- to understand the photosphere and dynamics of RSGs where about half of the atomic elements are produced ; 

- to explain the (unknown) mass-loss process which releases material to the surrounding medium and contributes 
heavily to the chemical enrichment of the Galaxy ; 

- to determine accurate abundances in distant galaxies from complex spectra by using proper model atmospheres. 

A previous attempt has been doneP to review the status of current or soon-to-be interferometers in terms of image 
synthesis. Interferometric data have been simulated from synthetic images in the K band (computed with the same image 
as in Fig. Q] and the same numerical codes used for this work). Reconstructed images of RSGs have been provided with 
cutting-edge code using the simulated data obtained with the VLTI facilities in a 6-8 telescopes configuration and with the 
earth rotation to provide a large (u,v) coverage. 

The reconstructed image of Betelgeuse shows in that case medium sized granules, similar to the original image smoo- 
thed to the interferometer resolution. There is, however, much more information in the original image (Fig.[TJ. Key features 
like the large central granule are very hard to see and the narrow inter-granular lanes are not seen at this angular resolution. 

Moreover, it appears that there is not apparent correlation of the intensity contrast between the smoothed and the 
reconstructed images. Then, the limitations of current image-reconstruction techniques lie in their ability to match exact 
intensities and to detect highly contrasted features. Indeed, the dynamic range is a function of the number of (u, v) points 
and of the accuracy of the visibility measurements.^ 

The timescale of the convective pattern is also crucial to constrain the simulations of red supergiant stars. A follow-up 
study at different epochs is probably needed at maximum 1 month apart in the optical wavelength where things change 
faster on timescales of days or even hours. This requires snapshot images and/or as many simultaneous measurements as 
possible in a short periodP^l 

It has been concluded^! that the complex structure of Betelgeuse requires a rich (u,v) plane sampling with a higher 
resolution and therefore the combination of the highest number of telescopes and of larger baselines. A good accuracy 
of the interferometric measurements is also mandatory by using for instance spatial filtering features and efficient fringe 
tracking systems. 
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Figure 2. Direct images of Betelgeuse with an annular array of 100 telescopes showing the crowding effect, which depends on the size of 
the clean field of view (CLF) compared to the object diameter (£>») : D, = 1 xCLF, D, = 1.5 x CLF, D, = 1.8 xCLF and D, = 2xCLF 
(from left to right). 
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Figure 3. Pearson correlation factor as a function of the ratio (star diameter)/(Clean FOV) to illustrate the crowding effect. 

2.4 Correlation criterium between images 

In this paper, we want to compare different images obtained with different array configurations and for different condi- 
tions of observations, so as to estimate the technical specifications of a hypertelescope to preserve an image of quality, 
regarding the astrophysical objectives of Betelgeuse. In a first approch, we have chosen as criterium the Pearson correla- 
tion coefficient,^ which allows to compare each image with the original one (Fig. [TJ. 

Pearson — — — 

where N p is the number of points in each image, X, (resp. Y[) is the intensity of the i™ point of the first (resp. second) 
image and X is the average of the N p intensities X/. This criterium equals to 1 in the best case (same images) and decreases 
down to in the worst case (no correlation between the two images). In a forthcoming paper, we want to define a more 
relevant criterium to estimate the restitution of both spatial and amplitude intensity of the structures in the image. 



3. IMAGING PROPERTIES IN THE CASE OF BETELGEUSE 



3.1 The aliasing effect and the crowding effect 

The spatial frequency coverage of a given array of telescopes is directly related to its imaging properties in the direct 
image plan. The highest spatial frequency or the maximum baseline B max corresponds to the smallest element of reso- 
lution in the image, called the Resel, with Resel - A/B max . The lowest spatial frequency or the minimum baseline B„,„, 
superimposes the lowest resolution, i.e. the largest imaging field of view, called the CLean Field of view : CLF - A/B„,j„, 
where an image can be formed properly. We can then estimate the number of Resels in the Clean FOV, which equals to 
Nr = {CLF /Resel) 2 = {B, nax /B m j„) 2 . In other words, the width of the Clean FOV equals to CLF = B max /B,„j„ expressed in 
number of Resels. 

In addition, the diameter of each telescope D, defines the portion of the sky seen by the interferometer, called the 
Coupled Field of view : CF = A/Di. A source located inside the CF (but outside the CLF) will be "coupled" to the 
interferometer by forming a halo of side-lobes in the detector plane. These sid e-lobe s induce photometric perturbations 
locally distributed in the image. This phenomenon is known as the aliasing effectP^H 

A consequence of the aliasing effect is the crowding effect. The crowding limit is defined as the maximum number of 
point-like sources Ns which can be contained inside the CLF : Ns < (B max /B m i„) 2 as defined above. Consequently, if the 
number of point-like sources inside the Coupled FOV (and inside the CLF) is larger than the number of Resels inside the 
CLF, the image inside the CLF will be saturated due to all the aliased sources. In the case of an extended source such as 
Betelgeuse, the crowding appears systematically if the star diameter exceeds the CLF width. Indeed, as the CLF is already 
completely filled by the central part of the star, the edge of the star outside the CLF will be aliased and will destroy the 
image inside the CLF. 

Figure|2]shows the direct images obtained with an annular array of 100 telescopes, which provides the same resolution 
in all directions on the sky, thanks to its circular symmetry. This figure illustrates the fact that it does not make sense 
to image a target that is substantially larger than the Clean FOV. This is confirmed by the figure [3] which shows that 
the Pearson correlation factor falls down above a certain value. Note that this value is not equal to 1, firstly because 1 
corresponds to the width of the image whereas the diameter of Betelgeuse is quite smaller, and secondly because the edge 
of the star is less bright than the center. Anyway, the crowding effect induces some smooth biases which increase when the 
CLF decreases. 

Thereafter, to avoid any problem of crowding, the CLF equals to the width of the image for all the following simula- 
tions in order to respect the sine-qua-none condition : D* < CLF. 

3.2 Influence of the array configuration on the interference function 

In a general sense for interferometry, the object on the sky is filtered spatially by the array of telescopes, so that the 
geometry of the array pattern and the number of sub-apertures determine the high angular resolution information sampled 
on the sky. That is why the choice of the array configuration is of a great importance, and depends on the scientific target. 

In the case of a hypertelescope, the image is biased by the halo of the interference functionfSl which induces a contrast 
loss in the image. The quality of the interference function of the array is simply related to the actual shape of the entrance 
pupil. If the input sub-pupils are distributed regularly, the densified pupil is almost complete, and the halo inside the Clean 
FOV reproduces the diffraction pattern of a large monolithic telescope covering all the sub-pupils. If the densified pupil 
shows gaps, additional diffraction figures are added to first ones. Thus, the halo is minimized by maximizing the densified 
pupil filling rate r„, with a regular pattern of the sub-apertures in the entrance pupilP-^ 

Figure|4]shows the case of 4 different geometries of arrays having the same number of telescopes (Nt ~ 100 here) and 
the same Clean FOV, i.e. the same minimum baseline. It appears clearly that for a given Clean FOV and a given number of 
telescopes, the OVLA configuration overcomes the other configurations, since we can see not only the large scale granules 
but also the inter-granular lanes. Indeed, the spatial distribution of the structures is barely disturbed, but the photometric 
distribution is biased by the halo of the PSF. 1 1 ELS A is clearly the worst one due to its non centro-symetric geometry. 
The images of CARLINA and KEOPS provide a better contrast between the bright and dark spots, but only the largest 
structures appear, which is mainly due to the fact that the maximum baseline is 3 times lower than for OVLA so as to keep 







Figure 4. Direct images of Betelgeuse (bottom) as a function of different array configurations of 100 telescopes having the same 
Clean FOV of 50 mas. Here is shown the output densified pupil of each array (top, from left to right) : annular configuration (OVLA), 
concentric multi-rings configuration (KEOPS), regular squared configuration (CARLINA), "Y" configuration (ELSA). For practical 
reasons, we use here the names in italic which corresponds to projects of large interferometric arrayi^even if they are already obsolete. 



the same minimum baseline. Note that as long as the object remains inside the Clean FOV and with a significant number 
of telescopes, the redundancy of the array has no influence (CARLINA or KEOPS provides the same image). 

In an other hand, the figure [5] shows the case of the OVLA and KEOPS configurations having the same Clean FOV 
and the same resolution, i.e. the same minimum and maximum baselines. This is made possible by adjusting the number 
of telescopes. It appears clearly now that KEOPS is the best one, thanks to its complete densified pupil which provides the 
best PSF quality with the lowest halo level and the best image with the highest contrast. '- However, KEOPS needs a huge 
amount of telescopes compared to OVLA. 

Figure [6] shows the Pearson correlation factor as a function of the array configuration (geometry and number of teles- 
copes). There is no significant gain of the image quality over 100 telescopes with OVLA. For KEOPS, the image quality 
increases with the number of telescopes along with the ultimate resolution. Firstly, both large and small scale structures are 
revealed. Secondly, the dynamic range improves so that the contrast between bright and dark structures tends to the same 
contrast in the original image (see Fig.Q}. 



3.3 Influence of the recombination mode on the diffraction envelope 

All the previous considerations (such as the resolution, the field of view, the crowding effect) are only related to the 
array, whatever the recombination scheme. We now consider the aspects related to the densified pupil combiner, for which 
different optical schemes are conceivable, 26 such as a fibered concept.^The standard scheme implies that the image should 
be multiplied by the diffraction envelope of an output densified sub-pupil. The width of the diffra ction envelope depends 
on the densification factor y, which is the scaling ratio between the output and the input sub-pupilsP^6'This factor can be 
adjusted between 1 < y < y max , where y = 1 is the Fizeau mode and y max = B m j n /Dj is the maximum densification when 
two output densified sub-pupils become tangent. We define also the Direct Imaging Field of view (DIF = A/((y - 1 ).£),-))) 
where the photons are concentrated inside the envelope. Note that most of the light falls within the Clean FOV in the 
maximum densification case (DIF « CLF)^ 




Figure 5. Direct images of Betelgeuse with the OVLA configuration (top) and the KEOPS configuration (bottom) having the same 
Clean FOV and an increasing number of telescopes. OVLA20 and KEOPS40 (left) have the same resolution (as 10 mas @ 745 nm). In 
the same way, OVLA36 and KEOPS 133 (middle left), OVLA63 and KEOPS408 (middle right), OVLA90 and KEOPS833 (right) have 
the same resolution (respectively « 5 mas, « 3 mas, ~ 2 mas @ 745 nm). 
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Figure 6. Pearson correlation factor as a function of the number of telescopes (left with OVLA and right with KEOPS). 

Due to the diffraction envelope, the intensification of the signal decreases from the axis as a Bessel function for circular 
apertures (or a pseudo-gaussian function in the fibered case^J), so that the intensity at the edge of the Clean FOV compared 
to the intensity on-axis is 40% lower in the maximal densification case and 10% lower with y = y„ MA /2p3Thus, a partial 
densification restitutes a more uniform photometry on the Clean FOV. In this paper, we neglecte this effect by considering 
a low densified image or a Fizeau image, since this effect is due to the combiner and is not linked to the array design. 



3.4 Object-Image-Hypertelescope relationship 

Figure [7] summarizes how to define the general design of the hypertelescope as a function of the scientific require- 
ments P3 
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Figure 7. Relationship between the astrophysical requirements, the direct image properties and the hypertelescope specifications. 
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Table 1. Technical parameters of a hypertelescope to image Betelgeuse in a field of view of 50 mas for different size of images expressed 
in number of Resels. 



The main dimensions of the object are the external diameter and the smallest Resel of interest. The maximum size of 
the object should not exceed the diameter of the clean field, which leads to the value of the minimum baseline of the array 
(CLF = A/B m j n ). The smallest Resel corresponds to the required resolving power, which imposes the largest baseline of 
the array (Resel = A/B max ). 

The complexity of the object determines the required number of Resels in the image, which depends on the number of 
telescopes and on the array geometry. The other aspect to be considered is the spectral density of the object (i.e. its spatial 
frequencies distribution) in comparison with the (u,v) plan coverage of the array on the range [B max /A, B mm /A], which 
provides a set of resolution samples A/B-, in the image for each baseline B,. 

The image quality, as well as the dynamic range, is directly related to the array geometry, which can be optimized 
by maximizing the densified pupil filling rate r Q . The limiting magnitude is directly related to the performances of the 
cophasing device, allowing long exposures. It depends also on the densification factor in the presence of read-out noise. 

To avoid aliasing and crowding, the object should not exceed the Clean FOV. Moreover, the diameter of the sub- 
apertures can be enlarged to avoid the aliasing of surrounding stars. In the case of Betelgeuse, the closest bright star 
(BD+07 1055B with a magnitude V 14.5 remote of 42 as) is so far that the aliasing cannot occur (as long as the telescope 
diameter is higher than 7 mm in the visible). Moreever, the presence of the circumstellar envelope around Betelgeuse can 
be problematic, but its luminosity remains negligible in the visible compared to the photosphere P2l 

Table[T]gives the main parameters of the hypertelescope to image Betelgeuse in a field of view of 50 mas. For instance, 
an image of 50x50 Resels with 1 mas of resolution at 745 nm of wavelength requires an array with B mm = 3.1 m and 
Bmax =154 m, 190 telescopes for OVLA and 3750 telescopes for KEOPS. The maximum densification factor is chosen 
to 10 so that the telescope diameter equals to D, = 0.31 m. The coupled FOV is 10 times larger than the Clean FOV 
(CF = 500 mas). A half-densification of y — y max /2 — 5 leads to a Direct Imaging FOV equal to DIF = 125 mas, 
where the photons are concentrated. Note that the ELT (diameter of 42 m) equipped with an aperture masking and a pupil 
densifier can provide a direct image of 13x13 Resels in the visible (@ 745 nm) and an image of 5x5 Resels in the infrared 
(@ 2.2 nm). 



4. TECHNICAL SPECIFICATIONS VERSUS CONDITIONS OF OBSERVATION 

4.1 Effect of the atmospheric pistons and cophasing aspects 

Figures[8]and|9]show the effect of the residual optical path difference (OPD) between the beams. As already knownpO 
the correction of the atmospheric piston is the hard point. An OPD standard deviation better than OPD > A/24 provides the 
best image quality and a cophasing specification of OPD > /1/8 is at least required. Otherwise, the small scale structures 
first and the large scale structures next are destroyed. This specification has also been demonstrated on a fibered testbedP^ 
Note that the integration of few hundreds frames allows to recover partially the spatial structures in the image if the OPD 
standard deviation remains higher than OPD > A/4. 

4.2 Effect of the photometric fluctuations 

Figures[lO]and[TT|show that the effect of the photometric fluctuations between the beams is not drastic for the hyperte- 
lescope. A differential photometry up to 50% does not really affect the image quality. A differential photometry up to 90% 
will generate a halo in the image but will preserve the structures in the image. This specification has also been demonstrated 
on a fibered testbed.^ 

4.3 Effect of the photon noise, the read-out noise and the background noise 

Figures[T2]and[T3]show the effect of different kind of noise in the image. A read-out noise is not critical as long as the 
signal to noise ratio remains higher than 10. The image becomes sensitive to the photon noise with a signal lower than 10 
photons per pixel. The image seems to be more sensitive to a background noise with a poisson signal lower than 100 pixel 
per pixel, but the thermal background and the zodiacal background remain negligible in the visible wavelength. 





Figure 9. Pearson correlation factor as a function of the standard deviation of the optical path difference (OPD) between the beams. The 
plots correspond to the average value (asterisk + dashed line) surrounded by the minimum and the maximum value (dotted line) of a 
data set of 100 frames, together with the value of the integrated data set (solid line). 



5. PLAN FOR IMAGE RECONSTRUCTION 

A more quantitative analysis of these images will necessarily rely on a more evolved process, involving image recons- 
truction techniques and taking into account the a priori knowledge of the PSF and/or the multi-chromatic information of 
the direct images recovered at different wavelength. 

A pseudo-deconvolution technique can be thought of by using the approximated formulation of the direct image as 
the convolution of the object by the interference function multiplied by the diffraction envelope.^ A first idea is that the 
diffraction envelope is well known (especially in the fibered case where the envelope is fixed and imposed by the fiber) and 
can be extracted as a flat-field, but the edge of the image is degraded compared to the center part, depending on the signal 
to noise ratio. However, the problem is that the strict convolution relationship is lost. The image and the PSF are in fact 
partially truncated, which is a problem for the classical methods of deconvolution. To overcome this problem, an hybrid 
method, based on likelihood maximization, reconstructing simultaneously the object and the PSF has been proposed.^ 

An other way is to use the Fourier information as currently done in interferometry for aperture synthesis imaging. Kno- 
wing that the hypertelescope concept is a generalization to N sub-apertures of the classical Michelson stellar interferometer 
with 2 telescopes, it has already been shown that for incoming monochromatic light, the only effect of pupil rearrangement 
is a shift of the high spatial frequency components of the object spectrum and that this shift can be easily corrected!^ Note 
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Figure 11. Pearson correlation factor as a function of the photometric fluctuations between the beams expressed in pourcents of 
maximum flux loss. Same legend as Fig. [9] 



that the spatial frequency peaks will be superimposed in the maximum densification case so that a partial densification 
(y < Jmaxl^-) and a non-redundant array are mandatory to preserve the Fourier information. Moreover, image reconstruc- 
tion can benefit from algorithms modelling the image to be reconstructed as a sparse combination of vectors in orthogonal 
representation bases like wavelets. A proposed approach treats the image reconstruction as a sparse approximation problem 
to include efficient and flexible prior geometric information for the reconstruction.^ 



6. CONCLUSION 

The astrophysical objectives and the array configuration are strongly dependent to design a hypertelescope. For a given 
number of telescopes, there is a trade-off between the resolution, the useful field of view and the contrast in the imageP-^ 
In a same way, for a given field of view and a given resolution (chosen as a function of the size of the object), there is 
a trade-off between the array geometry and the number of required telescopes to optimize either the (u,v) coverage (to 
recover the intensity distribution) or the dynamic range (to recover the intensity contrast). 

To obtain direct images of Betelgeuse with a hypertelescope, a regular layout of telescopes (like KEOPS) is the best 
array configuration to recover the intensity contrast and the distribution of the granulation cells but it requires a large 
number of telescopes (several hundreds or thousands). An annular configuration (like OVLA) allows a reasonable number 
of telescopes (lower than one hundred) to recover the spatial structures but it provides a low-contrast image. The large scale 
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Figure 13. Pearson correlation as a function of the read-out noise (left), the photon noise (middle) and the background noise (right). 
Same legend as Fig. [9] 



structure remains accessible with a limited number of telescopes (few tens) but the small granulation cells needs in any 
case many telescopes. 

To compensate the huge number of telescopes required in our simulations, a hypertelescope with few apertures can 
benefit from the earth rotational] and from the ability to re-configure telescope locations, although spotted rotated stars 
generally change in time too quickly to make this last solution efficient. Spectral dispersion can also be helpful to increase 
the (u, v) coverage by wavelength synthesis. 23 These aspects should be explore in future studies. 

The snapshot images provided by a hypertelescope are highly suitable to follow the evolution of the short timescale 
convective motions. The snapshot images allow also to feed the entrance plane of an integral field spectrometer to combine 
low, medium and/or high spectral resolution with the spatial resolution. But a relevant evaluation of the performance of the 
hypertelescope to study the complex astrophysical phenomena of the red supergiants requires anyway the development of 
a post-processing of the densified image. 

Finally, the photometric fluctuations are not critical (A photometry < 50%) contrary to the residual piston requirements 
(OPD < /1/8) which requires the development of an efficient cophasing system to fully exploit the imaging capability of a 
hypertelecope. 
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